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Abstract 

Muon bremsstrahlung photons converted in front of the DELPHI main tracker 
(TPC) in dimuon events at LEP1 were studied in two photon kinematic ranges: 
0.2 < Ery < 1 GeV and transverse momentum with respect to the parent 
muon p T < 40 MeV/c, and 1 < E 1 < 10 GeV and p T < 80 MeV/c . A good 
agreement of the observed photon rate with predictions from QED for the muon 
inner bremsstrahlung was found, contrary to the anomalous soft photon excess 
that has been observed recently in hadronic Z° decays. The obtained ratios 
of the observed signal to the predicted level of the muon bremsstrahlung are 
1.06 ± 0.12 ± 0.07 in the photon energy range 0.2 < E 7 < 1 GeV and 1.04 ± 
0.09 ± 0.12 in the photon energy range 1 < E 1 < 10 GeV. The bremsstrahlung 
dead cone is observed for the first time in the direct photon production at LEP. 
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1 Introduction 

Recent observation of anomalous soft photon production in hadronic Z° decays col- 
lected in the DELPHI experiment at LEP1 [1] has demonstrated the persistence of the 
soft photon anomaly found earlier in several fixed target experiments with high en- 
ergy hadronic beams, [2-6]. The photon kinematic range was defined in [1] as follows: 

0. 2 < Ej < 1 GeV, p T < 80 MeV/c, the p T being the photon transverse momentum 
with respect to the parent jet direction. Though the reaction e + e~ — > Z° — > hadrons 
presents a distinct mechanism of hadron production as compared to [2-6], the observed 
soft photon production characteristics were found in [1] to be very close to those reported 
in [2-6], both for the measured production rate and for the observed ratio of the rate to 
the inner hadronic bremsstrahlung. The latter was expected to be the main source of the 
direct soft photons in kinematic ranges under study (see [7-9]), while the observed sig- 
nals were found in [1-6] to be several times higher than the bremsstrahlung predictions. 
No theoretical explanation of this excess is available so far; reviews of the theoretical 
approaches to the problem can be found in [10,11] (see also the references [13-33] in [1]). 

From the experimental analysis, given a similarity of the soft photon production char- 
acteristics in both classes of experiments, the conclusion was drawn in [12] that the 
excess photons are created during the process of hadronization of quarks, i.e. their origin 
is strongly restricted to reactions of hadron production. If this ansatz is correct, a good 
agreement should be found between theory and experiment for the direct soft photon 
production in reactions of pure electroweak nature. What is the experimental situation 
in this field? 

The electron inner bremsstrahlung in e + e~ collisions (initial state radiation, ISR) was 
an important (and rather inconvenient) effect at LEP, with which all the LEP experiments 
had to contend. No deviation of the ISR characteristics from those expected from theory 
was observed, either at Z° or at high energy (see for example the DELPHI studies [13]). 
Therefore the situation with the electron inner bremsstrahlung can be considered as 
showing a nice agreement between theory and experiment. 

On the other hand, tests of QED with the muon inner bremsstrahlung which appears 
as final state radiation (FSR) in e + e~ — > fx + fi~ events were scarce at LEP. There were 
only two studies of photon production in Z° — > fi + fi~ events at LEP1 [14,15] and a single 
study of e + e~ — > events at LEP2 [16] 1 . All these studies aimed at the separation 

of rather hard photons, isolated from the neighbouring muon. So, the DELPHI analysis 
of final state radiation from muons at LEP1 [14] was restricted to the photon kinematic 
range of E 1 > 2 GeV, > 5°, i.e. to the transverse momenta with respect to the muon 
direction p? > 174 MeV/c. In [16] the minimum value of the angle 9^ was increased 
up to 15° (keeping the same photon energy threshold), tripling the minimum photon p T . 
The OPAL analysis at LEP1 [15] used photons of £ 7 > 0.9 GeV and 9 m > 200 mrad, 

1. e. the photon transverse momenta with respect to the muon direction were pr > 179 
MeV/c. Thus, an analysis of the muon inner bremsstrahlung in the soft photon kinematic 
range close to that analyzed in [1] is completely missing at LEP. This motivated us to 
study the reaction 

e + e _ — > Z° — > yU + /i~n7, n > 1 (1) 

at LEP1 in a photon kinematic range similar to the one analyzed in [1] (with the photon 
transverse momentum being defined now with respect to the parent muon direction). In 
addition to the low energy (LE) band of 0.2 < E 1 < 1 GeV explored in [1], a higher 

Outside the LEP experiments, a few studies of the muon inner bremsstrahlung have been done, see [17] and references 
therein. 
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energy (HE) band of 1 < E 1 < 10 GeV was also used in the analysis, being restricted 
however to the photons of small transverse momentum with respect to the parent muon 
direction, p T < 80 MeV/c. The pr range of the LE band chosen for the definition of the 
bremsstrahlung signal was taken narrower in this work as compared to that in [1] , namely 
Pt < 40 MeV/c. This choice was motivated by the fact that the photon angular variable 
used in this analysis, the photon polar angle relative to the parent muon direction, can 
be measured much more accurately as compared to the angular variable used in [1], the 
photon polar angle relative to the parent jet direction, and this confined most of the LE 
bremsstrahlung photons down to the mentioned Pt range. 

The results obtained in this study are presented both uncorrected and corrected for 
the photon detection efficiency The presentation of the uncorrected results is motivated 
by their better statistical accuracies and smaller systematic uncertainties in the absolute 
photon rates. 



2 Theoretical predictions for the muon inner 
bremsstrahlung 

In electroweak reactions like (1) the inner bremsstrahlung is a process of direct pho- 
ton production calculated via purely QED machinery. The production rates for the 
bremsstrahlung photons from colliding e + e~ (ISR) and from final [i + /i~ pairs (FSR) in 
the pt range under study can be calculated at once using a universal formula descending 
from Low [8] with a modification suggested by Haissinski [18]: 

where K and k denote photon four- and three-momenta, P are the 4-momenta of the beam 
e + , e~ and the muon involved, and p M is the 3-momentum of the muon; pi± = Pi — {n-pi)-n 
and n is the photon unit vector, n = k/k; r\ = 1 for the beam electron and for the outgoing 
fi + , f] = —1 for the beam positron and for the outgoing and the sum is extended 
over both beam particles and the parent muon (formula (2) is presented in the form of 
the photon production rate per muon); the last factor in the integrand is a differential 
production rate of the parent muon. 

As can be seen, formula (2) is of the lowest (leading) order in a. Higher order radiative 
corrections to it can be evaluated using exponentiated photon spectra in the LE and HE 
bands. In the accepted regions of low pr the effects of the exponentiation were found to 
be rather small, as considered in Sect. 6.2. 

To a great extent, formula (2) is used in this paper specifically to enable a compar- 
ison with the corresponding formula applied for the calculation of the inner hadronic 
bremsstrahlung in hadronic decays of Z° [1] (cf. the analogous formulae in [2-6]): 

^ - wte-J d PN %™*(p i K)(p j iod<K...<ppN (3) 

— * _ 

where K and k denote again photon four- and three-momenta, P are the 4-momenta of 
the beam e + , e~ and iV charged outgoing hadrons, and pi ... p^ are the 3-momenta of the 
hadrons; i] = 1 for the beam electron and for positive outgoing hadrons, rj = — 1 for the 
beam positron and negative outgoing hadrons, and the sum is extended over all the A^ + 2 
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charged particles involved; the last factor in the integrand is a differential hadron produc- 
tion rate (when calculating the photon production rate per jet only hadrons lying in the 
forward hemisphere of a given jet enter the sum). Calculations performed with formulae 
(2,3) show that the inner bremsstrahlung rate from one muon is approximately equal, in 
the kinematic region under study, to the predicted inner hadronic bremsstrahlung from 
a whole hadronic jet of a Z° hadronic decay. To a great extent, this is a consequence of 
the coherence of the photon radiation from the individual radiation sources, the charged 
hadrons produced in the fragmentation process. 

The contribution of the ISR to these rates is small, being below 1% in the photon 
kinematic range chosen for the analysis. This smallness is easy to understand: although 
the ISR from electron/positron beams is much more intense than the ISR from hadron 
beams in experiments [2-6], where it contributed a significant amount to the detected 
photon rate, all the extra photons in an experiment with colliding e + e~ are emitted at 
very small polar angles with respect to the beam directions, with the angular distribution 
peaking at 6 7 = Vs/T, where T is a beam Lorentz factor (r = 0.89 x 10 5 at the Z° peak), 
thus yielding few photons in the barrel region used in our analysis. 

The muon bremsstrahlung radiation (FSR) has the same angular behaviour of the 
photon production rate versus the photon polar angle relative to the parent muon direc- 
tion (the photon production angle, 6> 7 ), with T being in this case a muon Lorentz factor. 
For the muons from Z° decays at rest the T = 4.3 x 10 2 corresponds to the peak position 
at 4.0 mrad. Note that the position of the peak does not depend on the bremsstrahlung 
photon energy, since the dependences of the photon production rate on the photon energy 
and the photon production angle are factorized in formulae (2,3). The turnover of the 
muon bremsstrahlung angular distribution at the peak value and its vanishing at 6 y — > 
is termed the dead cone effect. This behaviour is illustrated by Fig. la where the initial 
part of the production angle distribution for the FSR of the reaction (1) is shown, gen- 
erated 2 with formula (2). The observation of the dead cone presents an experimental 
challenge requiring a highly accurate apparatus; the angular resolution of the opening 
angle between the measured muon and photon directions which is necessary for the ob- 
servation of the muon bremsstrahlung dead cone at LEP1, has to be of the order of 1—2 
mrad. 

3 Experimental technique 

3.1 The DELPHI detector 

The DELPHI detector is described in detail elsewhere [19,20]. The following is a brief 
description of the subdetector units relevant for this analysis. 

In the DELPHI reference frame the z axis is taken along the direction of the e~ beam. 
The angle 6 is the polar angle defined with respect to the z-axis, $ is the azimuthal 
angle about this axis and R is the distance from this axis. 

The TPC, the principal device used in this analysis, was the main tracker of the DEL- 
PHI detector; it covered the angular range from 20° to 160° in and extended from 
30 cm to 122 cm in R. It provided up to 16 space points for pattern recognition and 
ionization information extracted from 192 wires. The TPC together with other tracking 
devices (Vertex Detector, Inner Detector, Outer Detector and Forward Chambers) en- 
sured a very good angular accuracy of the muon track reconstruction, which is a part of 
the overall angular resolution for the photon production angle. The distribution of the 

2 The Monte Carlo data set of dimuon events described below was used as the input of the generation. 
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opening angles between the generated and reconstructed muon directions is shown in Fig. 
lb; it can be characterized by the distribution mean of 0.42 mrad and its r.m.s. width of 
0.37 mrad, which restricts 90% of the entries within the 0—1 mrad interval. 

The identification of muons was based on the muon chambers (MUC) surrounding 
the detector, the hadron calorimeter (HCAL) and the electromagnetic calorimeter (High 
density Projection Chamber, HPC), as described in [21]. 

The Monte Carlo (MC) data set used in this analysis was produced with the DYMU3 
generator [22]. Higher order radiative corrections to the reaction (1) total cross section 
were accounted for via the exponentiation procedure implemented in the generator. The 
generated dimuon events were passed through the DELPHI detector simulation program 
DELSIM [20]. 

3.2 Detection of photons 

Photon conversions in front of the main DELPHI tracker (TPC) were reconstructed 
by an algorithm that examined the tracks reconstructed in the TPC. A search was made 
along each TPC track for the point where the tangent of the track trajectory points di- 
rectly to the beam spot in the i?$ projection. Under the assumption that the opening 
angle of the electron-positron pair is zero, this point represented a possible photon con- 
version point at radius R. All tracks which have had a solution R that was more than 
one standard deviation away from the main vertex, as defined by the beam spot, were 
considered to be conversion candidates. If two oppositely charged conversion candidates 
were found with compatible conversion point parameters they were linked together to 
form the converted photon. The following selection criteria were imposed: 

• the <3> difference between the two conversion points should be at most 30 mrad; 

• a possible difference between the polar angles 9 of the two tracks should be at most 
15 mrad; 

• at least one of the tracks should have no associated hits in front of the reconstructed 
mean conversion radius. 

For the pairs fulfilling these criteria a x 2 was calculated from AO, A$ and the differ- 
ence of the reconstructed conversion radii AR in order to find the best combinations in 
cases where there were ambiguous associations. A constrained fit was then applied to 
the electron-positron pair candidate which forced a common conversion point with zero 
opening angle and collinearity between the momentum sum and the line from the beam 
spot to the conversion point. 

The photon detection efficiency, i.e. the conversion probability combined with the 
reconstruction efficiency, was determined with the hadronic MC data since the converted 
photon sample in dimuon events was insufficient statistically for such a determination. 
The efficiencies were tabulated against three variables: E 7 , 7 (the photon polar angle 
to the beam), and 0~,tk (the photon opening angle to the closest track). The efficiency 
varied with the energy from zero at the 0.2 GeV detection threshold up to 4 - 6% at 
E 7 > 1 GeV, depending on the two other variables (for details see [1]). 

In order to reduce a possible difference in the reconstruction of the converted pho- 
tons in the MC and the real data (originating from the bias in the detector material 
distributions in the two data sets and from a possible distinction in their pattern recog- 
nition results) the recalibration procedure described in [1] was implemented, with the 
recalibration coefficients obtained with hadronic events. 

The angular precision of the photon direction reconstruction was studied using the 
dimuon MC events and was found to be of a Breit-Wigner shape, as expected for the 
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superposition of many Gaussian distributions of varying width [23]. The full widths (Fs) 
of the A0 7 and A$ 7 distributions were 2.3 ±0.1 mrad and 1.9 ±0.1 mrad, respectively, 
for the combined 0.2—10 GeV interval (Figs, lc, Id). The full width of the distribution of 
the difference A# 7 between the generated and reconstructed muon-photon opening angles 
(which is the difference in the production angle # 7 defined in Sect. 2 and therefore repre- 
sents the overall angular resolution of the current analysis) was found to be 2.1 ±0.1 mrad 
(Fig. le), thus providing a possibility for the observation of the muon bremsstrahlung 
dead cone. Moreover, one can improve essentially this raw resolution, though at the price 
of a loss of 50% of the converted photon statistics, by requiring the photon energy to 
exceed 1 GeV and the conversion radius to be greater than 25 cm. With these tighter 
cuts 1.4 mrad resolution (the full width) was achieved and used in a particular case which 
required a high angular resolution and is described below (Sect. 7.3). 

The accuracy of the converted photon energy measurement was studied also with the 
dimuon MC events. In both energy bands it was at the level of ±1.5% (the Breit-Wigner 
full width about 3%); this is illustrated by Fig. If where the distribution of the relative 
difference between the generated and reconstructed photon energy is plotted for the LE 
photons. The resolution was checked with events of the (hadronic) real data by comparing 
the 7T° peak width of the 77 mass distribution from these data to the analogous one from 
the MC. 

4 Data selection 

4.1 Selection of dimuon events 

The data selection was done under standard cuts aimed at the separation of dimuon 
events (cf. [14,21]) which are described below. The consecutive application of these cuts 
reduced the MC sample of dimuon events by factors indicated in parentheses: 

• the number of charged particles N c h had to be within the interval of 2 < N c h < 5, 
and the two highest momentum particles had to have p > 15 GeV/c (0.894); 

• the polar angles of the two highest momentum particles had to be within the interval 
of2O°<0<16O° (0.962); 

• the impact parameters of the two highest momentum particles had to be less than 
0.2 and 4.5 cm in the and z projections, respectively (0.993); 

• no additional charged particles with momenta greater than 10 GeV/c were allowed, 
unless the fastest particle had a momentum greater than 40 GeV/c (0.999); 

• the acollinearity of the two highest momentum particles had to be less than 10° 
(0.989); 

• the two highest momentum particles had to be identified as muons using either the 
muon chambers (MUC), the hadron calorimeter (HCAL), or the electromagnetic 
calorimeter (HPC), by requiring associated hits in the muon chambers, or by energy 
deposition in the calorimeters consistent with a minimum ionizing particle (0.825). 

The total reduction factor for the MC events was 0.696. 

A total of 122 812 events of real data (RD) was selected under these cuts and compared 
to 373 918 selected MC events corresponding, after the normalization of the equivalent 
MC luminosity to the integrated RD luminosity, to 121 000 expected events. 
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4.2 Selection of photons 

The standard selection of converted photons was done under the following cuts: 

• only converted photons with both e + ,e~ arms reconstructed were considered; 

• 20° < 7 < 160°; 

• 5 cm < Rconv < 50 cm, where R con v means conversion radius; 

• 200 MeV < £ 7 < 10 GeV. 

384 and 1097 converted photons were found using these cuts in the real data in the 
LE and HE energy bands, respectively. Of these, 127 and 265 photons are in the selected 
p T regions: p T < 40 MeV/c for the LE band and p T < 80 MeV/c for the HE band. 

For a particular analysis done to scrutinize the dead cone effect (described in Sect. 
7.3), the photon energy was required to be between 1 and 10 GeV, and the conversion 
radius to be between 25 and 50 cm. 

5 Backgrounds 

The following background sources within the fi + fi~ event sample were considered: 

• External muon bremsstrahlung: 

the bremsstrahlung radiation from muons when they pass through the material of 
the experimental setup. 

• Secondary photons: 

when a high energy photon (of any origin) generates an e + e~ pair in the detector 
material in front of the TPC the pair particles may radiate (external) bremsstrahlung 
photons, which can enter our kinematic region. 

• "Degraded" photons: 

higher energy converted primary photons with degraded energy measurement due to 
the secondary emission of (external) bremsstrahlung by at least one of their electrons. 

DELSIM was invoked to reproduce these processes in the MC stream. 
Collection of background photons (all dubbed as External Brems) in the MC stream 
was done if any one of the following conditions was satisfied: 

• a given photon was absent at the event generator level, i.e. in the DYMU3 event 
record; 

• a given photon, found in the DYMU3 event record, migrated from outside a selected 
Pt region into that region due to the energy degradation. 

26.0 ± 2.9 and 61.5 ± 4.5 background photons (normalized to the RD statistics) were 
found in the selected pt regions: px < 40 MeV/c for the LE band and px < 80 MeV/c 
for the HE band, respectively. 

The background from Z° — > t + t~ events was estimated using the MC data produced 
with the KORALZ 4.0 generator [24] and passed through the full detector simulation and 
the analysis procedure. The r + r~ contamination of the dimuon event sample was found 
to be 1536 ± 20 events (1.3 % of the dimuon sample), which contain zero photons in the 
LE band and 1.3 ± 0.6 photons in the HE band, of which 0.3 photons would be in the 
Pt < 80 MeV/c region. In what follows, this background was neglected. 

The cosmic ray background was estimated from the real data, studying events which 
originated close to the interaction point, but outside the limits allowed for selected events. 
In both energy bands its contribution to the photon rates was below 0.1%. 
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The background from Z° — > e + e~ events tested with the BABAMC generator [25] 
with the full detector simulation was found to be vanishingly small. The same is valid 
for the 4-fermion backgrounds Z° — > e + e~fi + {i~ and Z° — > fi + fi~qq tested with events 
produced with generators [26,27]. 

6 Systematic errors 

6.1 Systematic uncertainties in the determination of the signal 

Since the converted photon sample in the dimuon event statistics collected by the 
DELPHI experiment during the LEP1 period was insufficient for the determination of 
the photon detection efficiencies and the recalibration coefficients, they were taken as 
being defined with hadronic events. Therefore it is worth to start the consideration 
of the systematic errors and their estimations with the uncertainties induced by these 
components of the analysis as they were determined in [1]. 

The uncertainty due to a difference in the photon propagation and conversion in the 
detector material in the RD and its simulation in the MC, and analogous difference in 
the pattern recognition, left after the recalibration procedure was applied (termed in [1] 
hardware systematics), was studied in [1] and evaluated to be 0.9% of the photon rate in 
the LE band and 2% in the HE band. 

The systematic error for the photon detection efficiency 3 , after the recalibration pro- 
cedure mentioned above being applied, is a purely instrumental effect originating from 
the choice of the binning of the variables used for the efficiency parametrization, resolu- 
tion effects, etc. In [1] it was found to range from 6% to 9% of the photon rate. These 
estimations were tested in the current study with the MC dimuon events by comparing 
the Pt spectra of the DYMU3 inner bremsstrahlung photons transported through the 
DELPHI detector by DELSIM with a subsequent simulation of their conversions, with 
the spectra of the same photons taken at the generator level and convoluted with the 
photon detection efficiency. In both energy bands the difference was below 5% which was 
the level of the test statistical accuracy. This means that the aforementioned error due 
to the detection efficiency is likely to be overestimated in [1], or it is really smaller in 
the muonic data, in particular, due to a better angular resolution and due to a narrower 
angular ranges in both energy bands. In what follows, the value of 5% is used as an 
estimate for the uncertainty of the efficiency calculation. 

The systematic errors originating from the influence of the pt resolution on the selec- 
tion cuts were estimated from runs with the reconstructed photon energy and production 
angle randomly shifted according to the appropriate resolution function (taking into ac- 
count the different angular resolutions in the LE and HE bands). The changes were found 
to be less than 0.3% of the photon rate in the LE band and less than 0.4% of the photon 
rate in the HE band. In what follows, the corresponding errors were neglected. 

The uncertainty of the background (BG) estimation is composed of the uncertainties 
coming from the DYMU3 generator, efficiency and hardware systematics, BG selection, 
and the procedure of the BG photon conversion simulation. The systematic error from 
the photon conversion simulation is considered to be equal to the systematic error of the 
photon detection in the MC stream, before the recalibration is applied, i.e. it can be 
approximated by the recalibration corrections, which were within 3-4%. The systematic 

li Notc that when dealing with the data uncorrected for the detection efficiency the efficiency error is relevant to the 
bremsstrahlung predictions only (since bremsstrahlung spectra have to be convoluted with the efficiencies in this case). On 
the contrary, when dealing with the corrected data the efficiency uncertainty has to be applied to the measured photon 
rates only. 
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errors due to efficiency and hardware in the background estimation have strong positive 
correlations with the analogous components of the systematic error in the calculation 
of the real data photon rates (indeed they are of the same relative amplitude, but the 
background errors have to be reduced by factors of 3.9 and 3.3 in the LE and HE bands, 
respectively, when entering the final systematic error, since the background rates within 
the corresponding p T intervals constitute 25.7% and 30.2% of the RD— BG photon rates 
in the corresponding energy bands). Ignoring, for the sake of clarity, these correlations we 
will consider all the background systematics components as independent and uncorrelated 
with the analogous components in the RD rates. Then, calculating the background 
systematic errors similarly to those for the RD and taking into account the suppression 
factors mentioned above, the systematic background uncertainties appear to be 1.4% and 
2.2% relative to the signal rate in the respective energy band in the case of the uncorrected 
data, and 1.9% and 2.7% in the case of the data corrected for efficiency. 
The above systematic errors are summarized in Table 1. 

Table 1. Systematic uncertainties (in % of the photon rates in the pt ranges below 
40 MeV/c and 80 MeV/c for the LE and HE photons, respectively) for the signal and 
the predicted muon inner bremsstrahlung. The total systematic error of each of the two 
energy band photon rates and signal-to-bremsstrahlung ratios is the quadratic sum of the 
corresponding individual errors, as quoted in Tables 2,3 below. 



Component 


Data uncorrected for 
the detection efficiency 


Data corrected for 
the detection efficiency 




LE band 


HE band 


LE band 


HE band 






Signal 






Recalibration 


0.9% 


2.0% 


0.9% 


2.0% 


Efficiency 






5.0% 


5.0% 


Background 


1.4% 


2.2% 


1.9% 


2.7% 




Predicted 


Bremsstrahlung 






Efficiency 


5.0% 


5.0% 






Formula (2) 


4.0% 


10.0% 


4.0% 


10.0% 



6.2 Estimation of the accuracy of the bremsstrahlung predic- 
tions 

The estimation of the accuracy of the bremsstrahlung predictions resulting from for- 
mula (2) was done by comparing FSR rates obtained with this formula and those delivered 
by the DYMU3 generator, in the corresponding p? ranges, as the difference between the 
predictions. In the LE band this uncertainty was about 4%, in the HE band about 10%. 
They are quoted in Table 1. 

These estimates agree well with the differences in the predictions for the muon inner 
bremsstrahlung rates obtained with formula (2), and those calculated with formulae which 
account for higher order radiative corrections, the calculations being performed with the 
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non-exponentiated photon spectrum [28] and with the exponentiated one [29,30]. In 
particular, the latter give 5.9% and 9.1% differences with formula (2) in the LE and HE 
bands, respectively. Note that when doing these calculations, parameter (5 which governs 
the bremsstrahlung photon spectrum [18] was obtained by integration of formula (1.2) in 
[29] applying our cuts, i.e. within rather narrow angular ranges varying as a function 
of the photon energy according to the p T cuts imposed in the corresponding energy band. 
The (3 values were found to be 0.0146 in the LE band and 0.0088 in the HE band, i.e. 
considerably smaller than (5 = 2a/n (lns/m^ — 1) = 0.0582, obtained by integration 
over all angles. The smallness of (3 reduces the difference between the bremsstrahlung 
predictions for the exponentiated and non-exponentiated photon spectra. 

7 Results 

Photon distributions for 6> 7 , pt and p\ are presented both for the data and the back- 
ground (left panels of Figs. 2-5), and for their difference (right panels of the figures). The 
latter distributions are accompanied by the calculated bremsstrahlung spectra according 
to Eq. (2) shown by triangles. 

To quantify the excess of the data over the background the difference between them, 
which represents the measured muon inner bremsstrahlung, was integrated in the px 
interval from to 40 MeV/c for the photons of the LE band, and from to 80 MeV/c 
for the photons of the HE band (these intervals correspond to the filled areas in panels 
d,f of Figs. 2-5), and the values obtained were defined as signals. However these pt cuts 
were not applied when filling the angular distributions displayed in Figs. 2-6 in order to 
keep these distributions unbiased. 

7.1 Energy band 0.2 < E 1 < 1 GeV, p T < 40 MeV/c 

Photon distributions, uncorrected and corrected for the photon detection efficiency, 
are displayed in Figs. 2 and 3, respectively. The results for the signal rate are given in 
Table 2 together with the predictions for the muon bremsstrahlung and their ratios. 



Table 2. The signal (RD— Background), the predicted muon inner bremsstrahlung 
(both in units of 10~ 3 7/yu) and their ratios in the px < 40 MeV/c range for the photons 
from the LE band. The first errors are statistical, the second ones are systematic. 



Value 


Data uncorrected for 


Data corrected for 




the detection efficiency 


the detection efficiency 


Signal 


0.412±0.048± 0.007 


25.9±4.0 ± 1.4 


Inner Bremsstrahlung 


0.388±0.001± 0.025 


23.30±0.01±0.93 


Signal/IB 


1.06±0.12±0.07 


1.11±0.17±0.07 



As can be seen from Table 2, the predicted and the measured muon bremsstrahlung 
rates agree well, within the measurement errors. The small differences in Signal/IB ratios 
between corrected and uncorrected data in Table 2 (and Table 3 below) arise from the 
non- uniformity of the efficiency reweighting factors. 
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7.2 Energy band 1 < E 1 < 10 GeV, p T < 80 MeV/c 

Photon distributions, uncorrected and corrected for the photon detection efficiency, 
are displayed in Figs. 4 and 5, respectively. The results for the signal rate are given in 
Table 3 together with the predictions for the muon bremsstrahlung and their ratios. 

As can be seen from Table 3, the predicted and the measured muon bremsstrahlung 
rates agree well, within the measurement errors. The smaller values of the corrected 
experimental and predicted bremsstrahlung rates in the HE energy band as compared 
to those in the LE band (while the energy range factor following from formula (2), 
\n(E™ ax / E™ m ) , works in favour of the HE band with an enhancement factor of 1.43) 
are explained by a higher attenuation of the rates induced by the px cut in the case of 
the bremsstrahlung photons from the HE band. 

Table 3. The signal (RD— Background), the predicted muon inner bremsstrahlung 
(both in units of 10~ 3 7//i) and their ratios in the < 80 MeV/c range for the photons 
from the HE band. The first errors are statistical, the second ones are systematic. 



Value 


Data uncorrected for 


Data corrected for 




the detection efficiency 


the detection efficiency 


Signal 


0.829±0.069± 0.025 


21.1±2.2± 1.3 


Inner Bremsstrahlung 


0.794±0.001± 0.089 


20.00±0.01 ±2.00 


Signal/IB 


1.04±0.09±0.12 


1.06±0.11±0.12 



7.3 Observation of the dead cone of the muon bremsstrahlung 

The distributions of the photon production angles with a fine binning (of 1 mrad 
bin width) are shown in Figs. 6a,b for the combined sample of the converted photons 
from both energy bands. The distribution obtained after background subtraction (Fig. 
6b) is accompanied by the calculated bremsstrahlung points. The displayed measured 
distributions are raw spectra, without any unfolding of the detector angular resolution; 
the bremsstrahlung spectra calculated with formula (2) were smeared instead by the 
resolution. We prefer to present the uncorrected measured distributions in order to 
demonstrate the independence of the obtained results on the correction procedure. 

As can be seen from the plots, the experimental points follow well the predicted 
bremsstrahlung distribution, showing a turnover at the expected bremsstrahlung peak 
position of 4 mrad. This is therefore an observation of the muon inner bremsstrahlung 
dead cone, for the first time in high energy physics experiments. The observation enriches 
the agreement between the experimental findings of the muon inner bremsstrahlung char- 
acteristics reported in this work and the QED predictions for the process. 

However a deeper insight into the bremsstrahlung pattern can be obtained when con- 
sidering, instead of the distribution dN 1 /d6 1 , the distribution dN^/dVL, where dfl is a 
solid angle element. Such a distribution is free of kinematic suppression at the polar an- 
gles 9 y approaching zero, and the remaining suppression of the photon production rate at 
very small angles is a purely dynamic effect, similar to that mentioned in Sect. 2 for the 
hadrons inside a jet, namely a destructive interference between the radiation sources, but 
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this time less straightforward, just between the muon "before" and "after" the photon 
emission 4 . The solid angle element dfl is proportional to <icos# 7 , which at small angles 
is, in turn, proportional to dO^. The position of the dN 1 /d6 2 distribution turnover is 
predicted to be at 6* = 1/T 2 (T = 430, see Sect. 2), i.e. at 9* = 5.4 x 10" 6 rad 2 . 

To observe this turnover, an improved angular resolution was required, achieved with 
the additional cuts (see Sect. 4.2) to be at the level of 1.4 mrad, as noticed in Sect. 3.2. 
The distribution dN^/dO 2 obtained with this resolution is shown in Fig. 6c, together with 
the bremsstrahlung predictions for this variable. Though the statistics are poor, the dip 
at 6* 2 < 5 x 1CT 6 rad 2 is visible in this distribution, revealing the dynamical dead cone of 
the muon inner bremsstrahlung. 

In order to estimate the statistical significance of this observation the following pro- 
cedure was undertaken. The initial part (about 20 bins) of the bremsstrahlung # 2 dis- 
tribution shown in Fig. 6c, with first two bins omitted, was fitted by a smooth curve 
(by a polynomial of 4th or 5th order). Then the fitting curve was extrapolated to zero, 
as shown in the figure giving the value of (5.64 ± 0.27) 7/5 x 10~ 6 rad 2 at the centre of 
the first bin of the distribution (the error reflects the variation in the fitting form and in 
the number of bins used in the fit). This value was assumed to represent the expected 
bremsstrahlung rate in the first bin of the distribution in a hypothetical situation when 
the bremsstrahlung dynamical dead cone is absent. The number of the real data photons 
in the first bin was 2 with the estimated background to be 0.66 ± 0.46, thus giving the 
signal value in this bin (1.34 ± 1.49) 7/5 x 10 -6 rad 2 . Assuming Poisson distribution 
for the signal photons these numbers correspond to the probability of the absence of the 
bremsstrahlung dead cone of less than 4%. 

8 Comparison with the hadronic soft photon analysis 

The main difference between the results of this analysis and the hadronic ones [1] is 
the absence of any essential excess of the soft photon production over the predicted inner 
bremsstrahlung rate reported in this study, contrary to the case for [1] where the observed 
soft photon rate was found to exceed the bremsstrahlung predictions by a factor of about 
4. The 95% CL upper limits on the excess factors which can be extracted from the results 
of this work are 1.29 in the LE band, and 1.28 in the HE band. 

Another distinction between the two analyses is an essential difference in the back- 
ground levels and in the possible systematic effects. However, the code transporting 
photons through the DELPHI detector and simulating their conversions in the detector 
material (DELSIM), the photon reconstruction algorithm and the determination of its 
efficiency, together with the recalibration procedure, were common to the two analyses. 
Thus the results of this work can be considered also as a cross-check of these procedures 
in the hadronic events study. On the other hand, the amount of dimuon events collected 
during the LEP1 period is considerably smaller than the number of collected hadronic 
events, due to a smaller Z° dimuon branching ratio (by a factor of 20). As a result, in 
the current analysis the statistical errors are either essentially higher than the systematic 
ones (in the LE band), or comparable to them (in the HE band), while in [1] the total 
uncertainties of the measured photon rates are dominated by systematic errors; neverthe- 
less it should be emphasized that the results of both analyses show clear signals of direct 
photons (even though the strength of the signal in [1] is not explained theoretically). 

4 In classical language, the radiation intensity into the solid angle dQ vanishes when three vectors: the muon velocity, 
its acceleration, and the radiation unit vector happen, in particular, to be parallel, see for example [31,32]. 
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9 Summary 

The results of the analysis of final state radiation in decays of Z° events at 

LEP1 are reported in this work. The radiation was studied in the region of small 
transverse momenta with respect to the parent muon, p? < 40 MeV/c in the photon 
energy range 0.2 < E 1 < 1 GeV (LE band), and px < 80 MeV/c in the photon en- 
ergy range 1 < E 1 < 10 GeV (HE band). The obtained photon rates uncorrected 
(corrected) for the photon detection efficiency were found to be, in units of 10~ 3 7//z, 
with the first error to be statistical and the second one systematic: a) in the LE 
band: 0.412±0.048 ± 0.007 (25.9±4.0 ± 1.4), while QED predictions for the muon in- 
ner bremsstrahlung were calculated to be 0.388±0.001 ± 0.025 (23.30±0.01 ± 0.93); b) in 
the HE band: 0.829±0.069±0.025 (21.1±2.2± 1.3), while the muon inner bremsstrahlung 
was calculated to be 0.794±0.001 ± 0.089 (20.00±0.01 ± 2.00). The obtained ratios of 
the observed signal to the predicted level of the muon inner bremsstrahlung are then 
1.06 ± 0.12 ± 0.07 in the LE band and 1.04 ± 0.09 ± 0.12 in the HE band (uncorrected 
rates are used for these ratios, as they possess a better statistical accuracy). Thus, the 
analysis shows a good agreement between the observed photon production rates and the 
QED predictions for the muon inner bremsstrahlung, both in differential (see Figs. 2-5) 
and integral (see Tables 2,3) forms. This is in contrast with the anomalous soft photon 
production in hadronic decays of Z° reported earlier in [1]. 

The bremsstrahlung dead cone is observed for the first time in the direct photon 
production in Z° decays in particular, and in the muon inner bremsstrahlung in the high 
energy physics experiments in general, also being in good agreement with the predicted 
bremsstrahlung behaviour. 
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Figure 1: a) Distribution of the final state radiation production angle in Z° — > 
events generated with formula (2); b) opening angle between the generated and recon- 
structed directions of a muon track; c) difference between generated and reconstructed 
photon polar angles 7 ; d) the same for the azimuthal angles $ 7 ; e) difference between 
generated and reconstructed fx'j opening angles 6> 7 , which illustrates the overall angular 
resolution of this analysis; f) difference between the generated and the reconstructed pho- 
ton energies in the photon energy range of 0.2 < E 1 < 1 GeV. The curves in Fig. lc-lf 
are the fits by Breit-Wigner forms (see text). 
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Figure 2: Photon distributions in the photon energy band 0.2 < E 1 < 1 GeV uncor- 
rected for the photon detection efficiency. The photon rates are given as the number of 
photons per 1000 muons per bin width of the distribution. Left panels: the data and 
background distributions for a) 7 , the photon production angle; c) photon pt] e) photon 
p\. Right panels, b), d), f): the difference between the data and the background for the 
same variables, respectively "Ext.Brems" corresponds to the background, "Int.Brems" 
corresponds to the muon inner bremsstrahlung predictions. The filled areas in panels d) 
and f) correspond to the signal integral (see text). The errors shown are statistical. 
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Figure 3: The same as in Fig. 2, corrected for the efficiency of photon detection. 
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Figure 4: Photon distributions in the photon energy band 1 < E 1 < 10 GeV uncorrected 
for the photon detection efficiency. The photon rates are given in number of photons 
per 1000 muons per bin width of the distribution. Left panels: the data and background 
distributions for a) 9 7 , the photon production angle; c) photon p T \ e) photon p\. Right 
panels, b), d), f): the difference between the data and the background for the same 
variables, respectively. The filled areas in panels d) and f) correspond to the signal 
integral (see text). The errors shown are statistical. 
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Figure 5: The same as in Fig. 4, corrected for the efficiency of photon detection. 
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Figure 6: Dead cone of the muon inner bremsstrahlung: a,b) as seen in the photon produc- 
tion angle distributions: a) the data and the background distributions; b) the difference 
between the data and the background; and c): the distribution of the photon production 
angle squared, obtained under tighter cuts that improve the angular resolution; the curve 
shows the fit of the bremsstrahlung distribution within 10~ 5 < 2 < 10~ 4 rad 2 by a 5th 
order polynomial extrapolated to the 1st bin of the distribution (see text). The errors 
shown are statistical. 



